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Recurrence plots exhibit line structures which represent typical behaviour of the 
investigated system. The local slope of these line structures is connected with a 
specific transformation of the time scales of different segments of the phase-space 
trajectory. This provides us a better understanding of the structures occuring in 
recurrence plots. The relationship between the time-scales and line structures are 
of practical importance in cross recurrence plots. Using this relationship within 
cross recurrence plots, the time-scales of differently sampled or time-transformed 
measurements can be adjusted. An application to geophysical measurements illus- 
. trates the capability of this method for the adjustment of time-scales in different 

^ ', measurements. 
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1 Introduction 



In the last decade of data analysis an impressive increase of the application 
of methods based on recurrence plots (RP) could be observed. Introduced 
by Eckmann et al. [1], RPs were firstly only a tool for the visualization of 
the behaviour of phase-space trajectories. The following development of a 
quantification of RPs by Zbilut and Webber [2,3] and later by Marwan et 
al. [4], has consolidated the method as a tool in nonlinear data analysis. With 
this quantification the RPs have become more and more popular within a 
growing group of scientists who use RPs and their quantification techniques 
for data analysis. Last developments have extended the RP to a bivariate 
and multivariate tool, as the cross recurrence plot (CRP) or the multivariate 
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recurrence plot (MRP) [5,6,7]. The main advantage of methods based on RPs 
is that they can also be applied to rather short and even nonstationary data. 

The initial purpose of RPs was the visual inspection of higher dimensional 
phase space trajectories. The view on RPs gives hints about the time evolu- 
tion of these trajectories. The RPs exhibit characteristic large scale and small 
scale patterns. Large scale patterns can be characterized as homogeneous, pe- 
riodic, drift and disrupted. They obtain the global behaviour of the system 
(noisy, periodic, auto-correlated etc.). The quantification of RPs and CRPs 
uses the small-scale structures which are contained in these plots. The most 
important ones are the diagonal and vertical/horizontal straight lines because 
they reveal typical dynamical features of the investigated system, such as range 
of predictability or properties of laminarity. However, under a closer view a 
large amount of bowed, continuous lines can also be found. The progression of 
such a line represents a specific relationship within the data. In this paper we 
present a theoretical background of this relationship and discuss a technique 
to infer the adjustment of time-scales of two different data series. Finally, an 
example from earth siences is given. 



2 Recurrence Plots 



A recurrence plot (RP) is a two-dimensional squared matrix with black and 
white dots and two time-axes, where each black dot at the coordinates {ti, ^2) 
represents a recurrence of the system's state x{ti) at time 



where m is the dimension of the system (degrees of freedom), £ is a small 
threshold distance, || ■ || a norm and ©(■) the Heaviside function. This def- 
inition of an RP is only one of several possibilities (an overview of recent 
variations of RPs can be found in [8]). 

Since R(ti,ti) = 1 by definition, the RP has a black main diagonal fine, the 
line of identity (LOI), with an angle of 7r/4. It has to be noted that a single 
recurrence point at (^1,^2) in an RP does not contain any information about 
the actual states at the times ti and t2 in phase space. However, it is possible to 
reconstruct dynamical properties of the data from the totality of all recurrence 
points [9]. 
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3 Line Structures in Recurrence Plots 

The visual inspection of RPs reveals (among other things) the following typ- 
ical small scale structures: single dots, diagonal lines as well as vertical and 
horizontal lines (the combination of vertical and horizontal lines plainly forms 
rectangular clusters of recurrence points). 

Single, isolated recurrence points can occur if states are rare, if they do not 
persist for any time, or if they fluctuate heavily. However, they are not a 
clear-cut indication of chance or noise (for example in maps). 

A diagonal line R(ti + T,t2 + r) = 1 (for r = 1 . . .1, where / is the length of 
the diagonal line in time units) occurs when a segment of the trajectory runs 
parallel to another segment, i.e. the trajectory visits the same region of the 
phase space at different times. The length of this diagonal line is determined 
by the duration of such a similar local evolution of the trajectory segments. 
The direction of these diagonal structures can differ. Diagonal lines parallel to 
the LOI (angle 7r/4) represent the parallel running of trajectories for the same 
time evolution. The diagonal structures perpendicular to the LOI represent 
the parallel running with contrary times (mirrored segments; this is often a 
hint of an inappropriate embedding if an embedding algorithm is used for 
the reconstruction of the phase-space). Since the definition of the Lyapunov 
exponent uses the time of the parallel running of trajectories, the relationship 
between the diagonal lines and the Lyapunov exponent is obvious (but this 
relationship is more complex than usually mentioned in literature, cf. [10]). 

A vertical (horizontal) line R(ti, ^2 + t) = 1 (for r = 1 . . .v, with v the length 
of the vertical line in time units) marks a time length in which a state does 
not change or changes very slowly. It seems, that the state is trapped for some 
time. This is a typical behaviour of laminar states [4] . 



4 Slope of the Line Structures 

In a more general sense the line structures in an RP exhibit locally the 
time relationship between the current trajectory segments. A line structure 
in an RP of length I corresponds to the closeness of the segment x{Ti{t)) 
to another segment x{T2{t)), where Ti(t) and T2{t) are two local time-scales 
(or transformations of an imaginary absolute time-scale t) which preserve 
that x(Ti(t)) ^ x{T2{t)) for some time t = 1...1. Under some assumptions 
(e. g. piecewise existence of an inverse of the transformation T{t), the two seg- 
ments visit the same area in the phase space) , a line in the RP can be simply 
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expressed by the time-transfer function 



(2) 



Especially, we find that the local slope b(t) of a line in an RP represents the 
local time derivative dt of the inverse second time-scale T2^{t) applied to the 
first time-scale Ti{t) 



This is the fundamental relation between the local slope b{t) of fine structures 
in an RP and the time scaling of the corresponding trajectory segments. Prom 
the slope b{t) of a line in an RP we can infere the relation 'd{t) between two 
segments of x(t) {'&{t) = Jb{t)dt). Note that the slope b(t) depends only on 
the transformation of the time-scale and is independent from the considered 
trajectory x{t). 

This feature is, e.g., used in the application of CRPs as a tool for the adjust- 
ment of time-scales of two data series [6,11] and will be discussed later. Next, 
we present the deforming of line structures in RPs due to different transfor- 
mations of the time-scale. 



5 Illustration Line Structures 

For illustration we consider some examples of time transformations for differ- 
ent one-dimensional trajectories f(t) (i. c. functions; no embedding). We study 
the recurrence behaviour between two segments /i and /2 of these trajectories, 
where we apply different time transformations to these segments (Tab. 1). In 
order to illustrate that the found relation (3) is independent from the under- 
lying trajectory, we will use at first the function f(t) — (Figs. lAl, Bl, 
CI etc.) and then f{t) = sin(7rt) (Figs. 1A2, B2, C2 etc.) as a trajectory. 
The local representation of RPs between these segments corresponds finally 
to cross recurrence plots (CRP) between two different trajectories/functions 
as will be mentioned later. 

Assuming that the second segment of a trajectory /2 is twice as fast as the 
first segment /i (Figs. lA), i.e., the time transformations are Ti{t) = t and 
T2{t) — 2t, we get a constant slope b — 0.5 by using Eq. (3). A fine in an RP 
which corresponds to these both segements follows '&{t) — 0.5 1 (Figs. lAl, 
A2). This result corresponds with the solution we had already discussed in 
[11] using another approach. In [11] we considered a simple case of two har- 
monic functions /i(t) = sin(Ti(t)) and f2{t) = sin(T2(t)) with different time 
transformation functions Ti = Lp ■ t + a and T2 = ip ■ t + (3. Using the inverse 
T^^ = and Eq. (3), we get the local slope of fines in the RP (or CRP) 



bit)^dtT2-'{T,{t))^d,m- 



(3) 
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Table 1 

Examplary time transformation functions Ti{t) and T2(i), the inverse T2^{t'), their 
corresponding slopes b{t) and time-transfer functions for lines in RPs shown in 
Fig. 1. 

Fig. Ti(t) T2{t) T^\t') h{t) m 

A t 2t 0.5 1' 0.5 0.5 1 

B t V0.2t' VoJt 

b = dtT2^ {T2{t)) — ^/if), which equals the ratio between the frequencies of 
the considered harmonic functions. 

In the second example we will transform the time-scale of the second segment 
with the square function T2{t) = bt^. Using Eq. (3) we get h{t) = and 

•&{t) — -\/0.2i, which corresponds with a bowed hne in the RP (Figs. IBl, B2). 
Since sin(7ra;) has some periods in the considered intervall, we get some more 
lines in the RP (Figs. 1B2). These lines underly the same relationship, but we 
have to take higher periodicities into account: ■i?(t) = \/0.2 /cvrt (/c e Z). 

The third example refers to a hyperbolic time transformation T2{t) — \l\ — t'^. 
The resulting line in the RP has the slope h{t) = — and follows ?9(t) = 
yl^^t^, which corresponds with a segment of a circle (Figs. ICl, C2). We 
can use this information in order to create a full circle in an RP. Let us 
consider a one-dimensional system, where the trajectory is simply the function 
f{T) — T{t), and with a section of a monotonical, linear increase T^^ — t and 
another (hyperbolic) section which follows Thyp = — \/r^ — t^. After these both 
sections we append the same but mirrored sections (Fig. 2A). Since the inverse 
of the hyperbolic section is T^^^ = ±\/r'^ — t"^, the line in the corresponding 
RP follows ^{t) = T^ypiTiinit)) = ±\/r^ — t^, which corresponds with a circle 
of radius r (Fig. 2B). 

An examplary data series from earth science reveals that such structures are 
not only restricted to artificial models. Let us consider the January solar inso- 
lation for the last 100 kyr on the latitude 44°N (Fig. 3A). The corresponding 
RP shows a circle (Fig. SB), similar as in Fig. 2B. From this geometric struc- 
ture we can infer that the insolation data contains a more-or-less symmetric 
sequence and that subsequent sequences are equal after a suitable time trans- 
formation which follows the relation T2'^{Ti) = \Jr'^ — t^. For instance, the 
subsequent sequences could be a linear increasing and a hyperbolic decreasing 
followed by a reverse of this sequence, a hyperbolic increasing and a linear 
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Fig. 1. Details of recurrence plots for trajectories f{t) whose sub-sections fi{t) 
and f2{t) undergo different transformations in time-scale (Tab. 1). Black areas 
correspond to times where fi{t) ~ /2(^)- The dash-dotted lines represent the 
time-transfer functions Note that these are not the entire RPs, only a small 
detail of them (an entire RP cannot contain only these structures - there are more 
features, like the line of identity (diagonal line from lower left to upper right) and 
a more or less symmetric plot around this line). 

decreasing part. Such bowed line structures can also be found in RPs applied 
to data from biology, ecology, economy. These deformations can obtain hints 
about the change of frequencies during the evolution of a process and may be 
of major interest especially in the analysis of sound data (an example of an 
RP of speech data containing pronounced bowed hues can be found in [12]). 

Whereas in the examples above only the second section of the trajectory un- 
dergoes a time transformation, in the last two examples (Figs. ID and E) 
the time-scale of the first section is also transformed. Nevertheless, the time- 
transfer function can be again determined with Eq. (2) too. 
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From these examples we can conclude that the line in a recurrence plot follows 
Eq. (2) and depends only on the transformations of the time-scale. 



A 



B 




Time t 



Fig. 2. Illustrative example of the relationship between the slope of lines in an 
RP and the local derivatives of the involved trajectory segments. Since the local 
derivative of the transformation of the time-scales of the linear and the hyperbolic 
sections (A) corresponds to the derivative of a circle line, a circle occurs in the RP 
(B). Recurrence plot is derived from the one-dimensional phase-space (no embedding 
used) . 

6 Cross Recurrence Plots 

The relationship between the local slope of line structures in RPs and the 
corresponding different segments of the same phase-space trajectory holds 
also for the structures in CRPs, 

CK{h,t2) = e(£ - \\x{h) - y{t2)\\). (4) 

which are based on two different phase-speice trajectories x{ti) and y{t2)- This 
relationship is more important for the line of identity (LOT) which then be- 
comes a line of synchronization (LOS) in a CRP [6,11]. 

We start with two identical trajectories, i. e. the CRP is the same as the RP of 
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Fig. 3. A corresponding structure found in experimental data: (A) the solar insola- 
tion on the latitude 44°N for the last 100 kyr (data from [13]) and its corresponding 
recurrence plot (B). The recurrence plot parameters were m = 1 and e = 2 (black) 
and e = 3.5 (gray). 

one trajectory and contains an LOI. If we now slightly modify the amplitudes 
of the second trajectory, the LOI will become somewhat disrupted. This offers 
a new approach to use CRPs as a tool to assess the similarity of two systems 
[5]. However, if we do not modify the amplitudes but stretch or compress the 
second trajectory slightly, the LOI will remain continuous but not as a straight 
line with an angle of 7r/4. The line of identity (LOI) now becomes the line 
of synchronization (LOS) and may eventually not have the angle n/A. This 
line can be rather bowed. Finally, a time shift between the trajectories causes 
a dislocation of the LOS, hence, the LOS may lie rather far from the main 
diagonal of the GRP. 

Now we deal with a situation which is typical in earth sciences and assume that 
two trajectories represents the same process but contain some transformations 
in their time-scales. The LOS in the CRP between the two trajectories can 
be described with the found relation (2). The function '&{t) is the transfer 
or rescaling function which allows to readjust the time-scale of the second 
trajectory to that of the first one in a non-parametrical way. This method 
is useful for all tasks where two time-series have to be adjusted to the same 
scale, as in dendrochronology or sedimentology [6]. 
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Next, wc apply this technique in order to re- adjust two geological profiles 
(sediment cores) from the Italian lake Lago di Mezzano [14]. The profiles cover 
approximately the same geological processes but have different time-scales due 
to variations in the sedimentation rates. The first profile (LMZC) has a length 
of about 5 m and the second one (LMZG) of about 3.5 m (Fig. 4). Prom both 
profiles a huge number of geophysical and chemical parameters were measured. 
Here we focus on the rock-magnetic measurements of the normalized remanent 
magnetization intensity (NRM) and the susceptibility n. 




Fig. 4. Rock-magnetic measurements of lake sediments with different time-scales. 
Corresponding sections are marked with different gray values. 

We use the time-series NRM and n as components for the phase-space vec- 
tor, resulting in a two-dimensional system. However, we apply an additional 
embedding using the time-delay method [15] (we do not ask about the phys- 
ical meaning here). A rather small embedding decreases the amount of line 
structures representing the progress with negative time [8]. Using embedding 
parameters dimension m = 3 and delay r = 5 (empirically found for these 
time-series), the final dimension of the reconstructed system is six. The cor- 
responding CRP reveals a partly disrupted, swollen and bowed LOS (Fig. 5). 
This LOS can be automatically resolved, e.g. by using the LOS-tracking al- 
gorithm as described in [11]. The application of this LOS as the time-transfer 
function to the profile LMZG re-adjusts its time-series to the same time-scale 
as LMZC (Fig. 6). This method offers a helpful tool for an automatic adjust- 
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ment of different geological profiles, which offers advantages compared to the 
rather subjective method of "wiggle matching" (adjustment by harmonizing 
maxima and minima by eye) used so far. 




Depth LMZC (m) 

Fig. 5. Cross recurrence plot between rock-magnetic data shown in Fig. 4. The 
dash-dotted Hne is the resolved LOS which can be used for re-adjustment the 
time-scales of both data sets. 



7 Conclusion 

Line structures in recurrence plots (RPs) and cross recurrence plots (CRPs) 
contain information about epochs of a similar evolution of segments of phase- 
space trajectories. Moreover the local slope of such line structures is directly 
related with the difference in the velocity the system changes at different 
times. We have demonstrated that the knowlege about this relationship al- 
lows a better understandig of even bowed structures occuring in RPs. This 
relationship can be used to analyse changes in the time domain of data series 
(e.g. frequencies), as it is of major interest, e.g., in the analysing of speech 
data. We have used this feature in a CRP based method for the adjustment 
of time-scales between different time-series. The potential of this technique is 
finally shown for experimental data from geology. 
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Fig. 6. Geological profiles after re- adjustment using the LOS which was found with 
the CRP shown in Fig. 5. Corresponding sections are marked with different gray 
values. 
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